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This paper develops models describing asexual and sexual replication in sporulating organisms. 
Replication via sporulation is the replication strategy for all multicellular life, and may even be ob- 
served in unicellular life (such as with budding yeast). We consider diploid populations replicating 
via one of two possible sporulation mechanisms: (1) Asexual sporulation, whereby adult organ- 
isms produce single-celled diploid spores that grow into adults themselves. (2) Sexual sporulation, 
whereby adult organisms produce single-celled diploid spores that divide into haploid gametes. The 
haploid gametes enter a haploid "pool", where they may recombine with other haploids to form a 
diploid spore that then grows into an adult. We consider a haploid fusion rate given by second- 
order reaction kinetics. We work with a simplified model where the diploid genome consists of only 
two chromosomes, each of which may be rendered defective with a single point mutation of the 
wild-type. We find that the asexual strategy is favored when the rate of spore production is high 
compared to the characteristic growth rate from a spore to a reproducing adult. Conversely, the 
sexual strategy is favored when the rate of spore production is low compared to the characteristic 
growth rate from a spore to a reproducing adult. As the characteristic growth time increases, or as 
the population density increases, the critical ratio of spore production rate to organism growth rate 
at which the asexual strategy overtakes the sexual one is pushed to higher values. Therefore, the 
results of this model suggest that, for complex multicellular organisms, sexual replication is favored 
at high population densities, and low growth and sporulation rates. 
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I. INTRODUCTION 

The emergence of sexual replication as the preferred 
replication strategy for complex multicellular organisms 
is one of the oldest and most important problems in evo- 
lutionary biology [1]. There are two broad theories for 
the selective advantage for sex [1-6] , both of which have 
more than one specific version. 

The first theory for the selective advantage for sex is 
the genetic repair theory, or, in the context of small pop- 
ulations, the Muller Ratchet theory [4, 6-10]. The ge- 
netic repair theory holds that sex evolved as a way to 
remove mutations from diploid genomes. By reproduc- 
ing via a haploid intermediate, defective copies of genes 
can be discarded, and functional copies of genes can be 
brought together via haploid fusion. 

In the context of finite populations, the genetic repair 
theory takes the form of the Muller Ratchet theory [7-9] . 
Briefly, Muller's Ratchet is a phenomenon whereby a fi- 
nite population will steadily accumulate mutations, and 
may eventually go extinct as a result. Sexual replication, 
by providing a mechanism to discard defective genes, can 
slow down or stop the Muller's Ratchet, and thereby pre- 
vent small populations from going extinct. 

The second theory for the selective advantage for sex 
is the adaptability theory [11]. This theory, which orig- 
inates with Weismann, has two versions: The Vicar of 
Bray hypothesis, and the Red Queen hypothesis. The 
Vicar of Bray hypothesis argues that sex increases vari- 
ability in small populations, making them more adapt- 
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able to changing circumstances. The name derives from 
an English cleric who was known for changing his religion 
according to circumstance [1]. Within the context of the 
idea that sex increases variation within a population, it 
is believed that host-parasite co-evolutionary dynamics 
drove the emergence of sex [6, 11]. This "arms race" ar- 
gument for the existence of sex is termed the Red Queen 
hypothesis. 

While the various theories for sex are not necessarily 
mutually contradictory, each of them is either incomplete 
or has difficulties. The Muller's Ratchet theory, for ex- 
ample, by relying on a small population size, suggests 
that sex should disappear in large populations. It is not 
immediately clear that this should be true, however, since 
many sexual organisms can attain seemingly fairly large 
population sizes. 

The Vicar of Bray and the Red Queen Hypotheses rely 
on a dynamic environment. However, there are numerous 
sexually replicating organisms that have remained essen- 
tially unevolved over millions of years, in what appear 
to be fairly stable environments. Therefore, it is not im- 
mediately clear that a dynamic environment is the main 
selective pressure driving the emergence of sex. 

The genetic repair theory is the theory with the broad- 
est acceptance among evolutionary biologists. That be- 
ing said, it could be argued that the theory is incomplete, 
because it is does not explain why sexual replication is 
disadvantageous in organisms such as bacteria [12]. Fur- 
thermore, it does not explain why, among the organisms 
that replicate sexually, some employ sexual replication 
merely as a stress response (such as Baker's yeast), and 
why others replicate exclusively via the sexual pathway. 

It should be noted here that the issue is not why re- 
combination is beneficial, since genetic recombination oc- 
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curs at all biological length scales [13]. The issue is 
why some organisms, such as bacteria, replicate asexually 
(and sometimes exchange genetic material with another), 
and why other organisms, such as humans, replicate via 
a haploid intermediate. 

In two recent papers [14, 15], Tannenbaum and Fonta- 
nari developed simple evolutionary dynamics models 
based on the sexual stress response in Baker's yeast. The 
first paper assumed a density-independent characteristic 
haploid fusion time, while the second assumed second- 
order haploid fusion kinetics. In both papers, the gen- 
eral result that emerged was that sex is favored when the 
characteristic haploid fusion time is small compared with 
the characteristic replication time. Therefore, the results 
of the models suggested that sex is favored in slowly repli- 
cating organisms, and that, all other factors being equal, 
sex is favored at high population densities. 

While the results of these models were broadly consis- 
tent with what is observed biologically, they are unsuit- 
able for considering additional features of sexual replica- 
tion, such as gamete differentiation into sperm and egg, 
and sex differentiation into male and female. Although 
many organisms replicate sexually, not all have distinct 
gametes, and, even among organisms that produce dis- 
tinct gametes (the anisogamous organisms), not all have 
distinct male and female sexes. There are therefore pre- 
sumably different regimes where the different types of 
sexual strategies will be advantageous. 

The previous two models, because they consider a sex- 
ual pathway whereby a given diploid cell splits into two 
haploids, do not reflect the fact that most sexually repli- 
cating organisms may continually produce gametes that 
may recombine to produce new organisms. The analo- 
gous process for asexually replicating organisms is sporu- 
lation. Therefore, before considering the regimes where 
asexual replication and various sexual replication strate- 
gies are advantageous, we first need to develop a formal- 
ism for describing the evolutionary dynamics of popu- 
lations which replicate by producing either single-celled 
spores or gametes. 

This paper is divided into four sections: In Section II, 
we develop a model describing asexual replication of or- 
ganisms that replicate by producing single-celled spores. 
In Section III, we consider the analogous process for sex- 
ual replication, where we assume that the gametes fuse 
via second-order kinetics. In Section IV, we numerically 
solve for the steady-state mean fitnesses of the various 
populations, and determine the various regimes where 
the different strategies are advantageous. We find that 
sexual replication is favored in organisms that produce 
offspring at lower rates. The critical offspring production 
rate increases as the characteristic organismal matura- 
tion time increases, or as the density of adult organisms 
increases. Therefore, the results of this paper appear to 
be broadly consistent with what is observed biologically. 
However, based on known scaling laws for organismal 
maturation times as a function of organism size, it is 
not clear that the results of our model imply that sexual 



replication will be favored in larger organisms. In Sec- 
tion V, we summarize our conclusions, and discuss future 
research directions. In particular, we discuss additional 
modeling that could resolve this possible inconsistency 
between our predictions for the preferred sexual strat- 
egy as a function of organism size, and what is observed 
biologically. 

We should note that the models we develop arc 
quasispecies-type models. Quasispecies theory has found 
application in a wide range of problems in evolutionary 
dynamics, including molecular, viral, and bacterial evo- 
lution, the immune response, and the emergence of can- 
cer [16-47]. The formalism is well-suited to studying the 
selective advantage for various types of replication strate- 
gies. 



II. ASEXUAL REPLICATION VIA 
SPORULATION 

In this section, we develop the evolutionary dynamics 
equations appropriate for describing asexual replication 
via sporulation. As discussed in the Introduction, while 
a model involving simple binary fission is appropriate for 
unicellular organisms, such an approach is inadequate for 
describing larger, multicellular organisms. The reason for 
this is that asexually replicating larger organisms begin 
their life cycle as a single cell. This cell develops into 
the adult organism, which then begins to produce single- 
celled spores at some rate (see Figure 1). 

Therefore, in order to properly study the regimes 
where asexual and various sexual replication strategies 
are advantageous, we first need to develop equations 
more appropriate for the replication dynamics of multi- 
celled organisms. 



A. Definitions 

The asexual replication of an organism (single or multi- 
celled) by sporulation occurs as follows: An immature 
organism grows to adult size. The adult then produces 
single-celled spores by budding. Each of these spores are 
immature organisms that then repeat the process and 
develop into adults on their own. 

For simplicity, we assume that the organisms have 
diploid genomes consisting of two chromosomes. The 
genome of each organism may then be denoted by 
{fi, er 2 }, where o\ and oi denote the respective base se- 
quences of each chromosome (we assume that each chro- 
mosome consists of a single strand of bases, so that a 
given chromosome o may be written as a — Sis 2 . . . s L . 
Here, L is the total sequence length, and each Sj is a base 
which is chosen from an alphabet of size S, where S 1 = 4 
for terrestrial life). 

Following standard practice in quasispecies theory, we 
make the simplifying assumption that there exists a sin- 
gle master sequence, denoted ctq, for which a given chro- 
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mosomc is functional. In this singlc-fitncss-pcak approx- 
imation, any chromosome a ^ ag is non-functional. A 
given genome may therefore be classified by the number 
of functional chromosomes it has, namely zero, one, or 
two. 

A genome with zero functional chromosomes is said to 
be of type (u,u), where u signifies that a chromosome is 
unviable. A genome with one functional chromosome is 
said to be of type (v,u), where v signifies that a chromo- 
some is viable. A genome with two functional chromo- 
somes is said to be of type (v, v). 

We assume that the growth of an immature organ- 
ism to adulthood is characterized by a first-order rate 
constant. We further assume that this rate constant is 
genome-dependent, since different organisms are differ- 
ently suited to the given environment, and so will reach 
maturity at different rates. We therefore let k vv , k vu , and 
k uu denote the first-order spore-to-adult rate constants 
for the (v,v), (v,u), and (u,u) organisms, respectively. 

We also assume that, once an organism reaches adult- 
hood, it produces spores at some fixed rate, which is 
again genome-dependent. The rates of per-organism 
spore production for the various genome types are de- 
noted w vv , lo vu , and lu uu . 

We let n a i.xy denote the number of asexual, imma- 
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ture organisms with genome (x,y), and n a 



denote 



the number of asexual, mature organisms with genome 
of type (x,y). 

During the sporulation process itself, the adult organ- 
ism replicates each of the chromosomes in the genome. 
It is assumed that the adult organism retains the origi- 
nal parent chromosomes, and only the two newly synthe- 
sized daughters segregate into the budding spore. This 
" immortal strand" mechanism is believed to occur in the 
adult stem cells of vertebrates and in budding yeast [48— 
50]. Since the spores are being produced by stem cells 
in the adult, it is reasonable to assume that a similar 
mechanism is at work here as well. 

Finally, the replication of each chromosome is not 
error-free. We let p denote the replication fidelity, de- 
fined as the probability that a v chromosome produces a 
v daughter. If we neglect backmutations, then a v chro- 
mosome produces a u daughter with probability 1 —p, and 
a u chromosome produces a u daughter with probability 
1. 



B. Mutation-selection equations 

The mutation-selection equations governing the evolu- 
tionary dynamics of the sporulating population are then 
given by, 
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n am is simply the total population of mature adults. We 
then define all population fractions with respect to this 
population number. Specifically, we define the popula- 
tion fractions x aq , rs — n aq , rs /n am , where q — m,i, and 
(r, s) is the genome type. Note that, in principle, some of 
the population fractions can be greater than 1, since they 
are not defined with respect to the total population, but 
rather with respect to the population of mature adults. 
Changing variables, we obtain, 
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where K a (t^ — k v vX a i VV -t- ftvuXai,vu ~t~ ^uuXai.uu is simply 
the mean fitness of the population. 



C. Steady-state mean fitness 

We may solve for the steady-state of this system of 
equations analytically. To begin, we first assume that 
K uu = 0, so that, the first three equations give, at steady- 
state, that, 

_ K a (t == oo) 

Xai.vv X arn ^ vv 

_ Rg(t = OO) 

Xai,vu X arrLjVU 

Xam,uu (3) 

so that, 

= [p 2 K vv LU vv -(K vv +K a (t = 00))K a {t = 00)]x am ^ vv (4) 

If x am ^ vv > 0, then we have, 

K a (t = OO) 2 + K vv R a (t = OO) - p 2 K vv LU vv = (5) 
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so that, 
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K a (t = 00) = K aA (t = 00) = \k vv [-1 + Jl + Ap 2 —] 

(6) 

If Xam.w = 0, then we have, 

= [pK vu LO vu {Kvu H~ K a (t — Oo))ft a (t = ^l^arn,™ (7) 

so, if x am . vu > 0, then we have, 



K a (t = 00) = K a fl(t = OO) = -K vu [-l + Jl + Ap — ] 

(8) 

Based on a stability analysis of the possible steady- 
states, it is possible to show that the actual value of 
R a (t = oo) is given by m&x{K a ,i(t — oo), K a ^(t = oo)}. 

Assuming that n vv > k vu , and that oj vv /k vv — 
w m /k to , we have that R(t — oo) — R a> i at p = 1. Defin- 
ing a — k vu /k vv and — lj vv /k vv , we obtain an error 
threshold transition at some p crit — p cr i t (a,/3), defined 
by the condition R a ,i(t = oo) = &a,2(t = oo). Below this 
replication fidelity, the effective growth rate of the (v, v) 
genomes is no longer competitive with that of the (v, u) 
genomes, and the result is the disappearance of the (v, v) 
genomes from the population. 



III. SEXUAL REPLICATION 

In this section we consider a sexual replication path- 
way, whereby a given diploid spore splits into two hap- 
loids, which may then fuse with other haploids in the 
population. 



For a haploid fusion rate governed by second-order re- 
action kinetics, the equations governing the evolutionary 
dynamics of the population are given by, 
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where V denotes the system volume in which the or- 
ganisms are present. In this model, we assume that the 
system volume changes in such a way as to maintain a 
constant density p of mature diploids (they are the fully 
grown organisms, so the total volume is dictated by the 
number of mature diploids present). We also define all 
population fractions with respect to the mature diploids, 
so that x s i,w = ^si,w / '^sm^vv j and x v — Ti v /n srn vv . 
In terms of the population fractions, we have, 
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where the mean fitness K s (t) — (1 /n sm , vv )dn sm ^ vv /dt = 
KwXsi,w Note that we do not need to include an equa- 
tion for x am>vv , since by definition x sm _ vv = 1. 



C. Steady-state mean fitness 



A. Definitions 

In the sexual replication model being considered here, 
a mature diploid with genome (v, v) produces smaller 
diploid spores at a rate given by oj vv . These diploid 
spores then divide into two haploid intermediates. The 
haploids of type v may then recombine with each other 
to form an immature diploid of type (v,v), which then 
grows to a mature diploid and begins the cycle again. 

In this model, we assume that only the v haploids may 
recombine with one another. Essentially, the u haploids 
are defective and cannot participate further in the repli- 
cation process. We assume that the haploid fusion rate 
is described by second-order kinetics characterized by a 
rate constant 7. We also assume that the haploids have 
a finite lifetime in the population, and decay with a first- 
order rate constant Kh- 

Because only viable haploids can recombine, the only 
diploids in the population are immature and mature 
diploids of type (v,v). We let n s i, vv , n sm , vv denote the 
number of immature and mature diploids of type (v,v), 
respectively. We also let n v denote the number of viable 
haploids in the population. 



Using the relationship between K s (t) and x si . vv , we ob- 
tain the steady-state equations, 

1 9 R s (t = 00) . . . 
= - 1P x 2 v - (1 + >-)R a (t = 00) 

= 2w vv p - jpxl - (kh + R s (t = oo))x v (11) 

The first equation may be solved for x v , giving, 
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Plugging into the second equation, and re-arranging, 
we obtain, 
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which may be solved to give, 



1 
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(15) 



Note that, at least in the limit of n vv /jp — > 0, we have 
K s (t = oo) > /c a (t = oo), with equality only occurring at 
p = 0, 1, or if = k vu . Therefore, in the model being 
considered here, sexual replication outcompctes asexual 
replication when the cost for sex (as measured by the ra- 
tio of the characteristic haploid fusion time to the char- 
acteristic growth time) is negligible. 



IV. COMPARISON OF THE VARIOUS 
REPLICATION MECHANISMS 

We may numerically solve for the mean fitness of the 
sexual replication mechanism, and compare the values 
of R a (t — oo)/k vv and R s (t = oo)/k vv . The population 
with the larger of the two values will be selected for in 
the given parameter regime. For simplicity, we assume 
that Hh = 0. In this case, the steady-state mean fitness 
of the sexual population is obtained by solving, 
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so that K s (t = qo)/k vv depends only on Kyy /7p and 
cu vv / k vv for a given value of p. 

The asexual mean fitness is given by R a (t — oo)/n vv = 
max{/t(t = oo) 0) i/K t)V , R(t = oo) 0)2 /k„ v }, where, 
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where we have assumed that lo vu /h vu — ui vv /k vv . 

For a given p we can determine, as a function of 
uj vv /k vv , the value of K vv /jp at which the sexual and 
asexual mean fitnesses are equal. This curve defines the 
boundary separating the regimes of dominance for the 
two strategies. 

We can analytically determine the behavior of this 
curve in both the low and high uj vv /k vv regimes. First, 
when oj vv /k vv is small, then we have, 
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and so, in any event, we have K a (t 

{u vv /k vv )xP, where \ = a ,P- 

Now, since we wish to determine the value of n vv /jp at 
which the sexual and asexual mean fitnesses are equal, we 



simply use the value of R a (t — oo)/k vv for R s (t = oo)/ k vv 
in Eq. (16). Since we are assuming lu vv /k vv is small, we 
can drop all higher-order terms in the numerator and 
denominator expressions, giving, 
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Now, when uj vv /k vv is large, then we have, 
R a ,i{t = oo) 



Kyy 

RaA 1 = °°) 
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and so, in any event, we have R a (t — oo)/k vv = 
P 1 / 2 x(w 1 , 1 ,/k 1 ,„) 1 / 2 , where \ = P 1 ^ 2 ,®. Following a sim- 
ilar procedure to the one used for small lo vv /h vv (only 
this time we drop the lowest-order terms since w vv /k vv 
is large), we obtain that the asexual and sexual mean 
fitnesses arc equal when, 
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We therefore find that, at low sporulation rates, the 
cost for sex, as measured by n vv /jp, must be made 
steadily larger as uj vv /k vv decreases in order for the asex- 
ual strategy to remain competitive with the sexual one. 
We find that this pattern indeed holds at all values of 
w V v/ K w The large oj vv /k vv behavior of this model also 
shows that, once the cost for sex drops below a critical 
value (given by 2(1 /\ 2 ~ 1) 2 )> then sexual replication is 
the dominant strategy for all values of u> vv /k vv . 

The cost for sex is measured by n vv /jp, because this 
quantity measures the ratio of the characteristic haploid 
fusion time (oc 1/jp) to the characteristic maturation 
time (oc 1/k vv ). What the results of this model sug- 
gest is that, as the sporulation rate decreases (as mea- 
sured by uj vv /k vv ), the critical cost for sex at which asex- 
ual replication overtakes sexual replication increases as 
well. If the cost for sex is sufficiently small, however, 
(< 2(l/x 2 — l) 2 ), then sexual replication outcompetes 
asexual replication for all sporulation rates. 

These results may be understood as follows: As the 
sporulation rate drops, the time interval between the 
production of diploid spores increases. Therefore, the 
additional fitness penalty incurred by having the diploid 
spores split and pay a time cost in finding other hap- 
loids with which to fuse decreases. The result is that the 
cost for sex has to be pushed to higher values before the 
strategy becomes disadvantageous. 

It is also interesting to note from the small ui vv /k vv 
expression for the critical value of n vv /^p that, as p de- 
creases, the critical cost for sex increases. This makes 
sense since, as the replication fidelity drops, the benefit 
of sexual recombination increases as well. 

Therefore, the results of this model suggest that, in 
the case of multicellular organisms that replicate via the 
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production of unicellular spores (or gametes, in the case 
of sexual organisms), sexual replication is favored in or- 
ganisms that sporulate slowly. This is essentially equiv- 
alent to the statement that sexual replication is favored 
in organisms that produce few offpsring. For such or- 
ganisms, the time cost associated with haploid fusion is 
comparatively small, and so it makes sense to replicate 
via a mechanism that provides the few offspring produced 
with the highest possible survival probability. 

Furthermore, as the maturation time of the organisms 
increases, or, as the population density increases, the sex- 
ual strategy becomes advantageous as well. Again, this 
makes sense, since both an increased maturation time 
and an increased population density reduce the time cost 
associated with sex, thereby leading to a selective advan- 
tage for the sexual replication strategy. 

V. CONCLUSIONS AND FUTURE RESEARCH 

This paper developed models describing asexual and 
sexual replication in sporulating organisms. Such mod- 
els are necessary for studying the selective advantages of 
asexual and sexual replication strategies in multicellular 
organisms. 

Under the assumption of second-order haploid fusion 
kinetics, we found that sexual replication is favored at 
low sporulation rates, or equivalently, at low offspring 
numbers, long organism maturation times, and high pop- 
ulation densities. These results make sense, since low 
offspring numbers, long organism maturation times, and 
high population densities all reduce the effective time cost 
associated with sex. In particular, for organisms that 
produce relatively few offpsring, it makes sense to invest 
a comparatively small cost in sex and thereby maximize 
the survival probability of each child. 

While the results of this paper appear to be consis- 
tent with actual behavior, a difficulty that arises is in 
correlating preferred replication strategy with organism 
size. It is well-known that the maturation time of an 
organism goes as TV 1 / 4 , where N is the number of cells 
in the adult organism [51, 52]. Therefore, as a function 
of N, we expect k vv to scale as A -1 / 4 . The population 
density of adults, however, should scale as 1/N, giving 
Kw/lP oc TV 3 / 4 . Therefore, the cost of sex goes to oo 
as N — > oo, and so, from the small uj vv /k vv criterion, 
we obtain that uj vv /k vv cx iV -3 / 4 for large N, so that 
u> vv cx 1/N. Therefore, as organism size increases, the 
rate of sporulation must scale as the reciprocal of the or- 
ganism size. This rate of decrease appears to be much 
too rapid to allow for the predominance of sex as the 



preferred replication strategy of larger organisms. 

To hopefully resolve this issue, future research will 
study the effects of mobile gametes on haploid fusion 
rates. In this vein, gamete differentiation (sperm/egg) 
and sex differentiation (male/female) will be studied as 
well. It is possible that, as organism size increases, these 
replication strategies reduce the time cost for sex to a suf- 
ficient amount that sex does indeed emerge as the pre- 
ferred replication strategy for larger organisms. While 
important work on gamete differentiation has been done 
by Duscnbery [53-55], a determination of the regimes 
where the various sexual replication strategies and asex- 
ual replication are advantageous has not yet been done. 

Another important factor that will need to be con- 
sidered is death. Thus far, our models do not assume 
that organisms eventually die. The neglect of this phe- 
nomenon could lead to an unrealistically large regime 
where asexual replication is dominant, though it is of 
course not yet clear if this is the case. 

The dependence of a and p on organism size and com- 
plexity will need to be considered as well. As organismal 
complexity grows, the size of the genome grows, and so it 
is likely that a and p decrease. The result is that u) vv /k vv 
must be pushed to higher values before asexual replica- 
tion becomes the dominant replication strategy, and this 
increase in the critical sporulation rate might be suffi- 
cient to result in sexual replication being the dominant 
replication strategy at larger organism sizes. 

Finally, a key assumption of our model is that the sex- 
ual organisms release gametes continuously. In reality, 
sexually replicating organisms generally store up gametes 
for a certain period, and then, during a mating season, 
collectively release these gametes into the surroundings 
(we are considering organisms that replicate in aqueous 
environments for this analysis). Thus, although the ga- 
mete production rate may be low, during these brief pe- 
riods when massive numbers of gametes are rapidly re- 
leased, the rate of haploid fusion is fast, thereby reducing 
the time cost associated with sex. This reduction in time 
cost could be quite large, which could significantly in- 
crease the size of the regime where the sexual strategy is 
dominant. In this case, a weaker dependence of lu vv on 
N may be necessary to ensure the selection for a sexual 
strategy, so that we do indeed obtain that sexual repli- 
cation is preferred as organism size grows. 
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